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Abstract 

The KATRIN experiment is designed to measure the absolute neutrino mass 
scale with a sensitivity of 200 meV at 90% C.L. by high resolution tritium 
/3-spectroscopy. A low background level of 10 mHz at the /3-decay endpoint is 
required in order to achieve the design sensitivity. In this paper we discuss a 
novel background source arising from magnetically trapped keV electrons in 
electrostatic retarding spectrometers. The main sources of these electrons are 
a-decays of the radon isotopes 219,220 Rn as well as /3-decays of tritium in the 
volume of the spectrometers. We characterize the expected background signal 
by extensive MC simulations and investigate the impact on the KATRIN 
neutrino mass sensitivity. From these results we refine design parameters for 
the spectrometer vacuum system and propose active background reduction 
methods to meet the stringent design limits for the overall background rate. 
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1. Introduction 

The Karlsruhe Tritium Neutrino (KATRIN) experiment is a next gener- 
ation, large-scale, tritium /3-decay experiment currently under construction 
at the Karlsruhe Institute of Technology (KIT); it will prospectively start 
taking data in 2015. KATRIN is designed to measure the effective electron 
anti-neutrino mass m J7e , defined as 
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where U e j denotes the Pontecorvo-Maki-Nakagawa-Sakata leptonic mixing 
matrix and m^ are the neutrino mass eigenstates JJ. The design sensitivity 
of KATRIN is 200 meV at 90% confidence level [2]. 

The experiment will use a model-independent technique based on the 
kinematics of tritium /3-decay. It will analyze the shape of the electron energy 
spectrum in a narrow region close to the tritium decay endpoint at E = 
18.6 keV. A non-zero neutrino mass reduces the maximum energy of the 
electron and changes the shape of the tritium /3-spectrum in the immediate 
vicinity of the endpoint. To reach the neutrino mass sensitivity, several 
criteria including high energy resolution, high signal count rates and low 
background must be fulfilled. 

In the 70 m long KATRIN setup (shown in figured]) a windowless gaseous 
tritium source (WGTS) of high stability and luminosity is combined with a 
large electrostatic retarding spectrometer of unsurpassed resolution j^. A 
magnetic guidance system adiabatically transports the electrons created in 
the tritium source towards the spectrometer where the energy analysis takes 
place. The spectrometer, working as an electrostatic filter, transmits only 
those electrons which have sufficient energy to overcome the retarding poten- 
tial. The transmitted electrons are then counted at a detector. By measuring 



the count rate for different filter voltages, the shape of the integrated energy 
spectrum can be determined. 

Since the spectrometer section must be essentially tritium-free, the tri- 
tium flow is reduced from the WGTS injection rate of 1.8 mbar • £/s down 
to a value of 10 -14 mbar • £/s at the end of the transport section. This un- 
precedentedly large suppression factor will be achieved by a combination of 
differential (DPS) and cryogenic pumping (CPS) elements 3 (]]. 

From the electron creation in the WGTS until the energy analysis in the 
central analyzing plane of the main spectrometer, the magnetic field drops by 
four orders of magnitude, collimating the electron momenta via the magnetic 
gradient force. This combination of Magnetic Adiabatic Collimation with 
Electrostatic filter, called the MAC-E filter principle, further described in 
section [21 allows for large solid angle acceptance, combined with high energy 



resolution 



7. 



In this paper we perform a detailed investigation of a novel background 
source arising from stored multi-keV electrons produced in /3-decays of tri- 
tium and secondary processes occurring during a-decays of the radon isotopes 
219 ' 220 Rn in the volume of the large main spectrometer. Due to the known 
magnetic bottle characteristics of a MAC-E filter for light charged particles, 
the electrons arising from nuclear decays inside the spectrometer volume are 
magnetically stored. With storage times of up to several hours, these parti- 
cles can significantly enhance the background level via ionization of residual 
gas. 

Despite the huge tritium retention factor, careful radio assays and use 
of low-activity components, the z/-mass measurements can be seriously dis- 
rupted by even single tritium /3-decays or a-decays of short-lived Rn-isotopes. 
This is due to exceedingly long storage times and the large number of back- 
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Figure 1: KATRIN experimental setup with main components a: rear section, b: 
windowless gaseous tritium source, c: differential pumping section, d: cryogenic pumping 
section, e: prc-spectromcter, f: main spectrometer, g: focal plane detector. Below, the 
magnetic field and the electric potential along the beam axis are displayed. In both 
spectrometers the MAC-E filter principle is applied: As the electric potential is increased 
to U rc t = —18.6 kV to filter the /3-electrons, the magnetic field drops from B c = 6 T 
to B c = 3 • 10~ 4 T, which collimates the electrons into a parallel beam with a flux of 
$ = 191 Tern 2 . 
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ground events resulting from one decay. We know of no other case in as- 
troparticle physics experiments where a single nuclear decay can continuously 
influence the measurements over a time period of several hours. 

Here we outline a detailed model of these processes validated by mea- 
surements at the much smaller pre-spectrometer 0, [loj| and use it to predict 
background rates and characteristics of the final KATRIN setup. Based on 
these results, we investigate implications on the neutrino mass sensitivity and 
demonstrate that the original KATRIN setup described in would result in 
background levels exceeding the required limits. To mitigate these problems 
we finally propose important design refinements and novel active background 
reduction methods. 

This paper is structured as follows: Section [2] briefly describes the fun- 
damental principles of the MAC-E filter as they pertain to the background 
described in this work. Section [3] gives a brief introduction to the software 
used to perform the MC simulations. In section HI the mechanisms of back- 
ground production will be outlined. In section [5j we discuss the expected 
event rates based on this background model for different vacuum scenarios 
in the spectrometer section, while the impact of this new background source 
on the KATRIN sensitivity will be discussed in section [6j 

2. MAC-E filter principle of KATRIN 

The WGTS produces electrons at a rate of ~ 10 11 Hz which are emitted 
isotropically and guided along magnetic field lines towards the spectrometer 
section, itself consisting of a smaller pre-spectrometer providing the option to 
filter out low-energy electrons, and a larger main spectrometer for precision 
filtering. The magnetic guidance of the electrons through the spectrometer 
is provided by a system of three superconducting coils (see figure [ID and an 
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external air coil system surrounding the main spectrometer. The retarding 
potential in both spectrometers is provided by inner electrodes constructed 



from wires, which allow for a 10 -6 precision of the filter potential 11] . The 
kinetic energy E^ [n of electrons entering the spectrometer section may be 
decomposed as 

Ekm = E± + E\\ , (2) 

where E± denotes the energy associated with the cyclotron motion and E\\ 
corresponds to longitudinal motion along a magnetic field line. Of E^ in only 
Eh is analyzed by the electrostatic filter. To achieve both high count rates and 
superior energy resolution, the initial E± component must be transformed 
into E\\ on the way to the central analyzing plane. This is achieved by 
the MAC-E filter principle, where, in the case of KATRIN, the magnetic 
field drops by four orders of magnitude from the entrance (or exit) of the 
spectrometer to its center (see figure [T]). By extending the reduction of 
the magnetic field strength over a length of about 10 m, a very smooth 
change of the magnetic field is assured, resulting in a fully adiabatic motion 
of the electrons. Due to this adiabaticity, the orbital magnetic moment /i of 
electrons is conserved. To first order, /x is given by 

E± 

a ~ ~ const. (3) 
\B\ 

The reduction of the magnetic field strength thus transforms the transversal 
energy at the edge (£j_) almost completely into parallel energy at the center 
of the spectrometer (Ef) 

I R | mm 

E% = E c kin -E c ± = E c kin - E e ± y— = E c kin - El ■ 5 • 10" 5 , (4) 

11 IT? max 

where quantities considered at the edge of the spectrometer carry superscript 
e, and those considered at center carry superscript c. It is E» which is 
analyzed by the electrostatic filter. 
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The MAC-E filter technique is as 
in direct neutrino mass experiments 



yet the most sensitive technique used 



l| , and a central design feature of the 



KATRIN experiment. The drawback, however, is that this magnetic field 
configuration inherently forms a magnetic bottle for light charged particles, 



since both ends of the spectrometer work as magnetic mirrors 
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3. Simulation tools 

The main principles of the MAC-E filter and its application in the KA- 
TRIN experiment can be understood analytically via the adiabatic approx- 
imation. However, in order to illuminate the role of the MAC-E filter as it 
applies to stored particle backgrounds and the complex, non-adiabatic situa- 
tions these entail, a precise and fast computational tool is required. The tasks 
of such a tool include the calculation of electromagnetic fields and particle 
trajectories to machine precision. This tool is provided by the simulation 
software Kassiopeia [jj], Q , which has been developed over the past years 



by the KATRIN collaboration. 
The trajectory calculations o 



Kutta methods described in 
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Kassiopeia are based on explicit Runge- 



19| . Electric and magnetic field calcula 



tions are performed via the zonal harmonic expansion 
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211 ] . In the case 



of electric fields, computations are carried out using the boundary element 



method [221 ] . Elastic, electronic excitation and ionization collisions of elec- 
trons with molecular hydrogen are included in the simulations; they are based 
on data and calculations in 23|-|3l|. The field, tracking and scattering simu- 



lations originate from FORTRAN and C codes developed between 2000 and 
2008 by one of us (F. C). 

In the framework of the investigations presented below, Kassiopeia was 
equipped with a selection of event generators including the /3-decay of tritium 
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and a-decays of different radon isotopes. The simulation o 
implemented using Fermi's theory of weak interactions 32 



tritium /3-decay is 



331 ] . Here we make 



use of a detailed tritium generator which includes the final state distribution 



of tritium 



35] and radiative corrections 36[ to the /3-spectrum, while 



shake-off effects at low energies are not yet included [371 ] . 

The modeling of electrons produced in radon a-decays includes processes 
described in more detail below, such as the creation of shake-off electrons pro- 
duced in the initial a-decay and conversion, shell reorganization and Auger 
electrons produced in the decay of the daughter polonium isotopes. The sim- 
ulation of these processes is based on data in the code Penelope 38J and the 



literature 39 



44]. 



The software has been validated by a number of associated measurements 



mostly performed as tes t experiments for KATRIN 



45j and cross-checked 



with other methods [4a, |47J as well as analytic calculations 48N53I] . The 
radon event generator in particular is validated through comparisons to pre- 



spectrometer 



0,0 



and independent measurements 



54]- 



4. Background production mechanism 

In this section, the basic ingredients of the background production mech- 
anism will be outlined. First, the nuclear decays of tritium and radon are 
examined in detail as sources of primary high-energy electrons in the keV 
range. Second, the electrons' dynamical behavior in a MAC-E filter and the 
mechanism of particle trapping and the conditions under which such stor- 
age may occur will be discussed. The final part of this section relates these 
processes to the observed background rates. 



S 



4-1. Nuclear decays as source of high- energy electrons 

The main source of keV-range primary electrons are nuclear decays. Of 
particular concern for the KATRIN experiment are tritium /3-decays and 
a-decays of the short-lived radon isotopes 219 ' 220 Rn. 

As a central design requirement of KATRIN Q , only an exceedingly small 
fraction of the order of 10 -14 of the tritium molecules injected into the WGTS 
will reach the spectrometer section. A small number of these molecules will 
decay there before being pumped out, thereby generating electrons with a 
continuous spectrum of up to about E . 

Electrons in a similar or higher energy range can be produced following 
nuclear a-decays (the primary a-particle as well as fluorescence X-rays are 
of no concern here). Due to the large pumping speed of the turbomolec- 
ular pumps (TMPs) connected to the main spectrometer volume and the 
correspondingly short pumping times (about 360 s), only short-lived radon 



isotopes are of importance here. There 



underground experiments like Gerda 55] . is not an issue for KATRIN due 



ore 222 Rn emanation, relevant for 



to its long lifetime {r-m = 5.51 d). On the other hand, the number of 219 Rn 
(T219 = 5.71 s) and 220 Rn (r 2 2o = 80.2 s) decays in the sensitive volume 
is not reduced significantly by pumping. Thus, these a-decays generate a 
background source which is distributed homogeneously over the entire spec- 
trometer volume (Vms — 1290 m 3 ). 

219 Rn arises from the 235 U actinide decay chain and emanates in small 
quantities primarily from the non-evaporable getter (NEG) material used 
for pumping the s pec trometers; for details of the Rn emanation from NEG 
material see 57 1. The stainless steel inner surfaces of the main and pre- 

spectrometer as well as auxiliary equipment such as ceramic insulators, glass 
windows, vacuum gauges and thermocouples emanate 220 Rn arising from the 
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Figure 2: Visualization of the different electron production mechanisms in 219 Rn 
a-decay. Details of the mechanisms are explained in the text. 

232 Th decay chain. The electron energy spectrum resulting from 219 Rn and 
220 Rn a-decays can be attributed to the processes of internal conversion, 
shake-off, shell reorganization and the Auger effect. 

Figure [2] is a sketch of the electron generation processes following radon 
a-decay. 219 Rn ( 220 Rn) decays into excited 215 Po* ( 216 Po*) states which then 
decay within a few picoseconds. If the wave function of a shell electron is 
non-vanishing at the nucleus, a conversion electron can be emitted in the 
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40j. This process 



de-excitation process, instead of radiating gammas 
is dominant for heavy nuclei, and in the case of polonium de-excitation, 
conversion electrons can reach energies of up to E e = 450 keV. 



The emitted a-particle can directly knock out shake-off electrons 
atomic shells. These electrons reach energies of up to E c = 80 keV 41 



rom the 
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44) 



Additionally, the emission of the a-particle results in a sudden, non-adiabatic 
change of the nuclear potential, which leads to the emission of predominantly 
two low-energy shell reorganization electrons from the outer shells which 
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Figure 3: Multiplicity of electrons following 220 Rn a-decays. The MC simulation 
with Kassiopeia based on the model described here is compared to the independent 
measurement in demonstrating the soundness of our background model. 



share an energy of about E c = 250 eV 43|. 

Subsequent to shake-off and conversion electron processes, which both 
may leave vacancies in the electron shell, Auger electrons can be emitted. The 



latter process often involves cascades of relaxations [38|, releasing multiple 
electrons with energies of up to E e = 20 keV. 

As an example of the complexity of the processes involved, figure [3] shows 
the electron multiplicity of 220 Rn a-decay as simulated for this investiga- 



tion and previously measured in an independent work , 



The simulation 



and measurements agree well, demonstrating the basic validity of our event 



generators. These generators are described in more detail in 



10|. 



4-2. The main spectrometer as magnetic mirror trap 

Due to the operating principle of the MAC-E filter, an electron produced 
in the center of the spectrometer is accelerated towards regions of low electric 
potential at the ends of the spectrometer, thereby moving from a region 
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of low to high magnetic field. Consequently, its longitudinal energy Eh is 
transformed into transversal energy E±. Depending on the starting angle 
and energy of the electron, the kinetic energy can be completely transformed 
into transversal energy so that the electron is magnetically trapped. 

Nevertheless, there are situations under which the storage conditions are 
broken: 

• Below a certain minimum transversal starting energy E± m , the electron 
cannot be magnetically trapped, since the acceleration by the electric 
field is too strong. In case of the main spectrometer, this minimal 
energy is Ef m = 0.93 eV. 

• Above a certain transversal starting energy E± ax , the electron's cy- 
clotron radius becomes larger than the radius of the main spectrome- 
ter, and therefore the electron directly hits the wall. For the reference 
field B m i n = 3 ■ 10 -4 T and the dimensions of the main spectrometer 
(0 = 10 m), this corresponds to E™ 3 * = 180 keV. 

• The motion of high-energy electrons in low magnetic fields can be non- 
adiabatic. Accordingly, the transformations of E± into E\\ and vice 
versa are no longer proportional to the change of the magnetic field, 
i.e. the angle of the momentum vector to the magnetic field line changes 
randomly. Therefore, non-adiabatic motion allows the electron to even- 
tually escape the magnetic mirror trap. 

4-3. Background production of stored high-energy electrons 

A primary high-energy (keV-range) electron which is trapped in the mag- 
netic bottle slowly cools down via ionization and electronic excitation colli- 
sions with residual gas molecules. Other energy loss mechanisms such as elas- 
tic scattering and emission of synchrotron radiation only play a minor role. 

12 




Figure 4: Calculation of the trajectory of a stored 3 keV electron in the main 
spectrometer with KASSIOPEIA. An electron produced via a nuclear decay can be 
trapped due to the magnetic mirror effect. Its motion is a superposition of cyclotron 
motion, axial motion and azimuthal magnetron drift. Here, the trajectory is shown for a 
time period about 100 /xs, whereas the actual storage time is much longer. 
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Due to the excellent ultra high vacuum (UHV) conditions (p = 10~ n mbar), 
collisions are rare, allowing a single electron to be stored for several hours. 
The hundreds of secondary electrons generated by ionizing collisions are 
mainly low-energy and typically leave the spectrometer on a rather short time 
scale of minutes. Accelerated by the retarding potential they hit the detector 
and thus produce a background in the narrow energy interval of the signal 
/3-decay electrons (the energy region-of-interest ROI is from 15 — 21 keV 2]). 

The total number of secondary electrons N s for a fixed primary energy 
given approximately by 

N s (E prim ) » (5) 

where u = 37 eV denotes the average energy of ion electron pair creation 
off H 2 for electrons in the keV-range 



5N 



and -Eprim represents the primary 



starting energy. For a realistic calculation of N s in our specific case, however, 
the following corrections need to be taken into account: 

• The high-energy secondary electrons themselves can be stored again 
and produce tertiary electrons, 

• at very high energies, electrons may leave the magnetic trap before 
being fully cooled down due to non-adiabatic effects, 

• stored electrons additionally lose energy by emitting synchrotron radi- 
ation. 

To incorporate these effects, we have carried out extensive simulations 
with the Kassiopeia simulation package, investigating the important pa- 
rameters of primary storage time t s (defined as the time between the cre- 
ation of the primary and the end of its trajectory) and number of secondary 
electrons N s as a function of E prim . The results are shown in figure [51 In 
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these simulations, eight different energies were selected in a range between 
10 eV - 100 keV. For each of these values, 10 3 electrons were started isotrop- 
ically in the main spectrometer. The simulation takes into account elastic 
scattering, excitation and ionization on H 2 at a pressure of p = lO^ 11 mbar, 
as well as non-adiabatic effects and synchrotron radiation. The average com- 
putation time of a single 10 keV stored electron and all its secondaries is 
^comp ~ 8 ■ 10 4 s on an Intel Xeon X5550 2.67 GHz processor. 

The trajectory of each electron (and all secondary electrons) was com- 
puted until it 

• leaves the spectrometer through the entrance or exit port, 

• or hits the spectrometer electrodes or vessel wall, 

• or was cooled down below the ionization threshold -Ethres = 15 eV. 

This cut-off parameter is motivated by the minor influence of electron ion- 
ization interactions below Ethrcs on our results. A detailed investigation of 
the processes below E thres will be described in a separate publication [59]. 



The results of our Kassiopeia simulations reveal a clear correlation of 
both background parameters with E pTira : a higher E py -i m implies a longer 
storage time (up to 10 h) and a higher multiplicity of secondary electrons. 
As a generic example, a 10 keV electron leads to the creation of ~ 300 
secondaries in a time period of 3 h, corresponding to a background rate of 
tb = 30 mHz in the energy ROI. For energies E prim < 30 keV, the average 
number of secondaries is a good means of estimation of the primary electron 
starting energy E prim . 

At energies above Eprim ~ 30 keV, the effects of non-adiabatic motion be- 
come more prominent. The large storage times and the considerable number 
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of secondary electrons underline the importance of the detailed investigations 
below. 



5. Expected background rates at the main spectrometer 

In this section, we use the background model described above to esti- 
mate the actual background rate to be expected at the final KATRIN setup. 
In doing so we examine four different vacuum strategies for the KATRIN 
spectrometer section to minimize the background rates. 

5. 1 . Sources of radon and tritium 

The pre-spectrometer measurements 9| initially revealed the dominant 
background source to be 219 Rn emanation from the 90 m of NEG strips (1.8 kg 
NEG material). Further sources of 219,220 Rn emanation were identified to be 
specific vacuum gauges and sensor instrumentation. After removal of the 
getter pump and the auxiliary components, a small number of events with 
radon-like characteristics were still observed, which we attribute to radon 
emanation from the inner surface of the stainless steel walls. 

The number of radon decays expected in the main spectrometer can be 
extrapolated from these measurements. In the main spectrometer a much 
larger NEG pump with 3000 m NEG strips (60 kg NEG material) will be 
used. The 219 Rn emanation for this batch was reduced by a factor of two 
through a special production process. However, since the decay series is not 
in secular equilibrium, the radon production rate increases slowly over time 



at a rate of 0.3 Bq/(kg • yr) [60 ] . 



As mentioned before, the remaining radon events after removal of the 
getter and auxiliary components are assumed to be caused by emanation from 
the walls. For the following discussion, we assume that this radon emanation 
rate scales with the respective spectrometer vessel surface Aj (surface area 

16 



■■■Fit 

a Simulation 



non-adiabaticity 



4? 



10 



10 2 



10 3 10 4 10 s 

Primary-electron start-energy [eV] 











MM 


Fit 


non-adiabaticity , 


i 




• Simulation 






1 1 1 Mill 




ii 










nil 










• * 







10 



10* 



10 J 10 4 10 b 

Primary-electron start-energy E . [eV] 



prim 



Figure 5: Average storage time t s (top) and average number of secondary elec- 
trons N s (bottom) as a function of the primary electron starting energy -Eprim. 

The error bars signify the standard deviation away from this average. The Kassiopeia 



simulations reveal a power law behavior scaling with t s ~ -E p rim ±0 ^ an< ^ N s ~ -E'prim 
for -Ep r i m < 30 keV. Above 30 keV the effects of non-adiabatic motion become dominant, 
as indicated by the shaded area, reducing the average storage time of the primary electron 
and the number of secondaries, respectively. 



(0.97±0.02) 
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of the pre-spectrometer: Apg « 25 m 2 , main spectrometer: Ams ~ 690 m 2 ). 
We assume a negligible emanation of radon isotopes from vacuum gauges and 
sensor instrumentation, as well as from the structural materials of the inner 



electrode system [11]. 



To calculate the number of tritium decays in the main spectrometer, we 
use the maximum allowed tritium flow from the WGTS into the pre-spectro- 
meter, approximately Qj? s ~* PS ~ 10 -14 mbar-^/s = 2.5 -10 5 molecules/s, as 
detailed in Q] . Additionally, due to the large number of adsorption/desorption 
processes in the transport section, we note that the gas flow will be a mixture 
of hydrated tritium (HT) and other hydrogen isotopologues, in the context 
of this work, however, we only consider T 2 . 

5.2. The vacuum system of the KATRIN spectrometer section 



611 ] is based 



The vacuum system of the KATRIN spectrometer section 
on two pumping strategies: TMPs to pump out noble gas atoms, such as 
radon, and secondly NEG pumps for pumping out hydrogen isotopologues, 
including tritium. Accordingly, for the NEG system, there is an inherent 
trade-off between increased tritium pumping capacity and enhanced radon 
emanation. 

Specifically, the vacuum system of the spectrometer section consists of six 



TMPs 62[ installed in pairs at the three pump ports at the detector-facing 
end of the main spectrometer and two smaller TMPs at the pre-spectrometer 
pump ports. The pump ports in the main spectrometer (pre-spectrometer) 
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are additionally equipped with 3000 m (180 m) of NEG strips 

To reduce the number of radon atoms reaching the main spectrometer 
volume from the NEG strips, LN 2 -cooled cryo-baffle systems were installed 
in front of each of the three pump ports. The excellent performance of this 
method to shield the sensitive flux tube of a spectrometer from the pump 
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port which houses the NEG strips has been demonstrated previously 65 1. 

The presence of the baffle, however, also results in a decrease of the hydro- 
gen (tritium) pumping speed. As will be shown in the following, an optimum 
solution is found for a configuration with 250 m additional getter strips in 
the pre-spectrometer and the installation of cryo-bafHes in all three pump 
ports of the main spectrometer. 

Table [1] summarizes the sources and corresponding reduction rates of 
radon atoms and tritium molecules for specific layouts of the KATRIN vac- 
uum system. 

5.3. Calculation of the expected decay rates 

The time dependent number of radon atoms N M (Rn) obeys the following 
differential equation: 

d7V MS (Rn) = _ _ n ms {Ru) _ N MS (Rn) gMS(Rn) + Q PS^M S(Rn) + ^M S(Rn) ) 

dt v v ' VMS s v ' s v ' 

total decay rate v inflow from PS emanation rate 

pump out rate 

(6) 

where AR n denotes the radioactive decay constant of the corresponding radon 
isotope, Vms stands for the volume of the main spectrometer and S' MS (Rn) 
is given by the sum of the available pumping systems for a particular radon 
isotope in the corresponding UHV scenario. The back-flow of radon from the 
main spectrometer to the pre-spectrometer can be neglected. In equilibrium 
we expect 

diV MS (Rn) = 

di 1 ' 

and one finds the number of radon isotopes iV MS (Rn) to be 

AT-dta) = [C— (Ha) + JS-fRa)] ■ A „„. yMs %. S (Rn) - < 8 » 

The number of tritium molecules iV (T2) in the main spectrometer is 
described similarly by 

7V MS (T 2 ) - NMS ^ . S MS (T 2 ) + Q PS ^ MS (T 2 ), (9) 



s v ' Vms 



total decay rate v inflow from PS 

pump out rate 



19 



Tabic 1: This tabic displays all important input parameters for the calculation of radon and 
tritium decay rates. We list the emanation rates of radon in the main spectrometer from 
3000 m getter strips with and without the baffle installed NEG , baf , NEG ) 

as well as from the wall (-E^aii)- As an °th er source, the inflow of radon and tritium 
from the pre-spectrometer for the cases of 180 m getter (Q^scT^neg) an( l 180 m plus 
additional 250 m getter strips (Qadd^250 m neg) are gi yen - Furthermore, the effective 
pumping speeds of the TMPs (S^f MPs ), and the NEG pump with baffle (S 3 M ^ m NEG + baf ) 
and without (S^^ m neg) are listed for tritium and radon. Finally, we list the radioactive 
decay constants (A). For all calculations a sticking coefficient of 0.8 for the cryo- baffle is 
assumed. The tritium inflow rates are to be understood as upper limits and are therefore 
given without an error estimation. The errors on radon emanation rates are propagated 
from measurements at the pre-spectrometer. 
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Tabic 2: Description of scenarios. The scenarios are differentiated in their selection of the 

amount of NEG strips in the pre-spectrometer and main spectrometer as well as the usage 

of LN2 cooled baffle. 

pre-spectrometer main spectrometer 

180 m NEG add. 250 m NEG 3000 m NEG LN 2 baffle 

Scenario 1 — — — — 

Scenario 2 / •/ / — 

Scenario 3 / / / / 

Scenario 4 / — / / 

where At 2 denotes the radioactive decay constant of T 2 and S ,MS (T 2 ) is 
given by the sum of the available pumping systems for T 2 for a specific UHV 
scenario. 

Using the input parameters summarized in table [IJ the decay rates are 
computed for four different scenarios outlined in table [2j 

• Scenario 1 completely avoids background from NEG correlated radon 
decay activity, 

• Scenario 2 primarily reduces background arising from tritium decay, 

• Scenario 3 optimally reduces both background rates arising from tri- 
tium and radon decays, see figure O, 

• Scenario 4 uses no additional getter in the pre-spectrometer. This sce- 
nario will be realized at the start-up of the spectrometer test measure- 
ments in 2012. 

Table [3] shows the nuclear decay rates corresponding to the different sce- 
narios (table EJ) expected in the main spectrometer. The table clearly demon- 
strates the importance of the NEG strips to reduce the number of tritium 



21 



Tabic 3: Expected number of nuclear decays in the main spectrometer for different UHV 

scenarios. The scenarios are described in detail in the main text. 

Activity [mBq] 

T2 219 RnNEG 219 R.nwaii 220 R.nwaii 

Scenario f 21.9 27.2 ± 27.2 67.9 ± 45.0 

Scenario 2 0.01 118.7 ± 29.5 27.2 ± 27.2 67.9 ± 45.0 

Scenario 3 0.02 0.12 ± 0.03 5.3 ± 5.3 1.4 ± 0.1 

Scenario 4 0.05 0.05 ± 0.01 5.3 ± 5.3 1.4 ± 0.1 



/3-decays in the main spectrometer volume as well as the non-negligible num- 
ber of radon decays even in the case of optimum passive shielding of the pump 
ports. With an overall decay rate of the order of a few mBq, primarily due 
to radon emanation from the inner spectrometer walls, a concise calculation 
of the resulting background rates is mandatory 

5.4- Expected background rates in different vacuum scenarios 

With the above information on the number of decays as well as the number 
of secondaries produced in each decay, the total expected background rate 
can be calculated. The average number of background events (Nb) in the 
energy ROI and in a time interval t longer than the storage time t > t s is 
given by 

(^) = ^ S -^ S E<^ MS X^ MS >, (10) 

i 

where i denotes the isotopes tritium and 219 ' 220 Rn, (N^) stands for the av- 
erage number of nuclear decays in a time interval t and (N ei ) represents the 
average number of electrons produced within one event. 

Two further factors in our background estimate have to be taken into 
account: the sensitive volume of the main spectrometer amounts to only 
6y S = 0.7 of the total volume, and only a fraction e^ s = 0.4 of all secondary 
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electrons produced in the main spectrometer will propagate towards the de- 
tector (the remaining 60% fly towards the source side). This asymmetry in 
the exit direction is due to the asymmetric magnetic field configuration (see 
figure [1]) . 

Assuming the partial pressures of tritium and radon to be constant over 
long time periods, the decay rate follows a Poisson distribution. The distri- 
bution of the number of secondary electrons is obtained by MC simulations 
of 10 3 tritium /3-decays and 219 ' 220 Rn a-decays each (based on the event gen- 
erators described in section [3]). 

On the basis of these considerations, we investigate in detail the influence 
of specific design modifications on the overall background rate. As an impor- 
tant example we briefly discuss the merits of additional getter strips in the 
pre-spectrometer. In figure El the trade-off between a reduced background 
from tritium and an increased radon-induced background due to additional 
NEG strips is clearly visible. 

For the tritium retention factor as stated in 2j, the tritium-induced back- 
ground is larger than the 219 Rn-induced one by about a factor of five if no 
additional getter is installed. Were the actual tritium retention factor to dif- 
fer from the reference value, figure |6] would allow adjustments to be made to 
the pre-spectrometer getter length. 

These two isotopes, however, are only part of the overall background 
picture, which is displayed in figure [7J for the four UHV scenarios listed in 
section 1531 From figure [7] it is evident that scenarios 1 and 2 result in overall 
background rates of ~ 1 Hz, thus exceeding the design criterion of 10 mHz 
by about 2 orders of magnitude. 

When comparing scenarios 3 and 4 one notes that the total background 
rate is almost identical. This is because the rate is largely dominated by 
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radon emanation from the inner surface of the main spectrometer, due to 
the excellent shielding of 219 Rn emanation from the pump ports by the LN 2 - 
cooled baffles. However, even when including these passive measures, the 
expected overall background rate of ~ 30 mHz still exceeds the design crite- 
rion by a factor of three, pointing to the need for additional active background 
reduction techniques. 

6. Impact of the background on the neutrino mass sensitivity of 
KATRIN 

As outlined above, radon emanation from the inner walls of the spectrom- 
eter and its structural materials may easily exceed the reference background 
level by a factor of three, and potentially, in case of larger than expected 
emanation rates, the background level would be correspondingly larger. 

In case of a Poisson-distributed background Nb g , the statistical uncer- 
tainty cr stat on the observable m| of KATRIN will scale roughly as ~ nK 6 



66(. However, the background arising from stored electrons is not of this type. 
Instead, the fluctuations of the rate are largely determined by the number of 
stored primary particles, which is small compared to the number of secon- 
daries reaching the detector. The count rate at the detector therefore shows 
rather large fluctuations which are not Poisson-distributed. Accordingly, the 
variance is determined by the variance of the number of stored primaries 
arising from nuclear decays. 

To investigate the impact of the background arising from stored electrons, 
a detailed model describing the background as a function of time over the 
full three years measurement time of KATRIN was implemented. The model 
is based on the full MC simulations described in section 15.41 and calcula- 
tions discussed in section 15.31 The statistical sensitivity m^ tat (90%C.L.) is 
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Figure 6: Estimated contribution of background rate originating from tritium 
and 219 Rn decays as a function of additional getter length in the pre-spec- 
trometer. As the plot shows, the radon contribution to the background increases with 
increasing amounts of getter material whereas the tritium contribution decreases. The 
optimum for a tritim retention factor as given in 2} is found at about 250 m additional 
getter in the prc-spectromcter volume. The rather shallow minimum in this case is due 
to the coincidence of almost identical rates of tritum and radon induced background. In 
case that the tritium inflow would be larger by an order of magnitude, more NEG strips 
would be required. The pre-spectrometer is able to hold up to 1000 m of getter. 
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Figure 7: Expected background rates for four different UHV scenarios described 
in table [2] and section 15.31 In this figure, resulting background rates are broken down 
into detail. For each scenario, the contribution from each distinct nuclear background 
source is displayed, illustrating the trade-offs and consequences inherent in implementing 
the various techniques. Scenarios one and two differ in their use of NEG to combat tritium 
background at the expense of introducing radon; the similarity of total rates here is purely 
coincidental. In comparison to scenarios one and two, three and four include a cryo-bafflc 
designed to trap radon emanating from NEG strips. 
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determined by fitting the theoretical integral /3-spectrum to 10 4 simulated 
KATRIN measurements. Each simulation assumes m Fe = eV and statisti- 
cally corresponds to three years of data taking. The width of the distribution 
of the fitted neutrino mass squared ml determines the statistical uncertainty 
cr s tat and thereby the neutrino mass sensitivity at 90% confidence level ac- 
cording to 



mf t (90%C.L.) = Vl.64-<r 8tat . (11) 

In a typical measurement schedule, the integral tritium /3-spectrum will be 
measured at 41 different retarding potentials. The overall measurement time 
at each potential is optimized to achieve the best neutrino mass sensitiv- 
ity Q for a background level of 10 mHz. During a measurement period of 
three years, a large number of scans of a few hours' duration through all 41 
measurement points will be performed. 

Figure [8] shows the statistical neutrino mass sensitivity at 90% confidence 
level as a function of the overall background rate (leaving all other contri- 
butions at their reference values). Here we compare a Poisson-distributed 
background (as used in 2j) to the background model as calculated in this 
work including nuclear decays, using a fixed scan time of t scan = 3 h and a 
pressure level of p = 10~ n mbar. 

When comparing both results, it becomes evident that statistical error 
increases significantly in case of a non-Poissonian background. It is this par- 
ticular feature of nuclear decays that necessitates the development of active 
background reduction techniques to realize the full physics potential of KA- 
TRIN. 

Moreover, our simulations revealed that m^ tat (90%C.L.) strongly depends 
on the actual time of a scan t scan and on the total pressure p in the main 
spectrometer. At scan times i scan much longer than the storage times t s , 



27 



the non-Poissonian nature of the background from nuclear decays becomes 
prominent. Since t s decreases inverse proportionally with p, an analogous 
effect is observed for higher pressures. Consequently, both larger values of 
t Bcan and larger values of p will result in a decrease of the neutrino mass 
sensitivity. Furthermore, in an ordered scanning mode, measurements at 
neighboring filter potentials are correlated due to the long storage times t S can- 
By scanning the 41 potentials in a random order, this correlation will be 
alleviated, and consequently the neutrino mass sensitivity can be improved. 
These interdependencies are visualized in figure [9] and [KB 

In summary, our investigations point to the following important facts: 

• The estimated neutrino mass sensitivity m^ tat has to take into account a 
detailed background model, the experimental scan mode and the UHV 
conditions of the spectrometers. 

• Backgrounds from nuclear decays feature large non-Poissonian rate fluc- 
tuations, which result in a decrease of the neutrino mass sensitivity 
m^ 

• The neutrino mass sensitivity m^ tat improves with better vacuum, with 
smaller scanning times, and with randomized scanning (instead of or- 
dered scanning). 

7. Conclusion and Outlook 

Due to their inherent electromagnetic design features, the KATRIN spec- 
trometers act as magnetic bottles for light charged particles. A primary 
electron in the multi-keV regime produced by a nuclear decay can thus be 
magnetically trapped over a time period of several hours during which it can 
produce several hundred secondary electrons. 
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Figure 8: Statistical neutrino mass sensitivity as a function of background rate 
for Poisson distributed background and for the background model including 
nuclear decays. For these simulations, a fixed scan time of i S can = 3 h and a pressure of 
p = 1 • 10~ u mbar was used. The dashed lines indicate the statistical sensitivity reached 
with a Poisson-distributcd background of 10 mHz (as stated in Q) and with the estimated 
background level of 30 mHz (see figure [7]) arising from nuclear decays of this work, if no 
active reduction methods are implemented. 
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Figure 9: Statistical neutrino mass sensitivity m^, tat as a function of pressure 
in the main spectrometer. For this simulation, a background level of 30 mHz and 
a scanning time of t sciin = 3 h was assumed. For a constant Poisson-distributcd back- 
ground (dashed line) the statistical sensitivity does not depend on the pressure. When 
the non-Poissonian background arising from nuclear decays is included, the experimental 
sensitivity m^, tat gets worse for higher pressures, as described in the text. The impact of 
the non-Poissonian nature of this background can be alleviated by scanning the potentials 
in random order (triangles), as compared to a fixed order (dots). 
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Figure 10: Statistical neutrino mass sensitivity m^, tat as a function of scan- 
ning time. For this simulation, a background level of 30 mHz and a pressure of 
p = 1 • 10~ n mbar was used. The presence of a non-Poissonian background results in 
a dependence of m^ tat on the scanning time in an analogous way as on the pressure, see 
figure [HI Again, when scanning the potentials in a random way (triangles) the neutrino 
mass sensitivity is improved, as compared to a scanning method with a fixed order (dots). 
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In this paper, we showed that nuclear decays of tritium migrating from the 
WGTS to the spectrometer as well as of 219 Rn and 220 Rn emanating from 
NEG material and structural components in the volume of the KATRIN 
main spectrometer can cause a background rate exceeding the design limit 
of 10 mHz. 

In an optimum scenario, using LN2-cooled baffles to shield the pump 
ports, as well as optimized combination of NEG strips and TMPs, a back- 
ground level of ~ 30 mHz is expected. 

Of major impact for the neutrino mass sensitivity m^ tat are the large 
rate fluctuations that this background exhibits. A statistical analysis with 
a detailed background model revealed a reduction of the statistical neutrino 
mass sensitivity from m^ tat = 0.16 eV at 90% C.L. (assuming a Poisson 
distributed background of 10 mHz) to m^, tat = 0.25 eV at 90% C.L. (assuming 
a realistic background model of 30 mHz). 

This result highlights the necessity for further developing; active back- 
ground reduction methods. In a separate publication 67| we describe the 



successful implementation of the electron cyclotron resonance (ECR) method, 
which offers great potential in reducing the background described here to a 
very low level. The upcoming measurements with the KATRIN main spec- 
trometer starting in the second half of 2012 will be of crucial importance to 
test this and other promising active background reduction methods. 
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